Microscopic model of a phononic refrigerator 
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We analyze a simple microscopic model to pump heat from a cold to a hot reservoir in a nanome- 
chanical system. The model consists of a one-dimensional chain of masses and springs coupled to a 
back gate through which a time-dependent perturbation is applied. The action of the gate creates a 
moving phononic barrier by locally pinning a mass. We solve the problem numerically using a non- 
equilibrium Green function technique. For low driving frequencies and for sharp traveling barriers, 
we show that this microscopic model realizes a phonon refrigerator. 



I. INTRODUCTION 

In atomic gases, techniques such as evaporative cooling 
can bring temperatures down to the submicrokelvin scale, 
allowing for the observation of quantum phenomena such 
as Bose-Einstein condensation. Progress in refrigerating 
condensed matter systems has been less spectacular but 
steady; a recent example is the use of active feedback for 
cooling nanomechanical cantilevered beams;-^ Other ex- 
perimental examples, most of them in electronic devices, 
have been reviewed in Refs. i^- On the theoretical side, 
mechanisms for quantum refrigeration by means of elec- 
tronic pumping with ac fields operating at low^ and high^ 
frequencies have been recently proposed. There are also 
several proposals for cooling nanomechanical systems by 
absorbing phonons with electrons i^ 

The challenge of quantum refrigeration is to keep the 
heat extracted from the cold sample higher than the 
energy dissipated into it by operating the cooling ma- 
chine. In electronic ac pumping, this condition can 
be achieved within the so-called "adiabatic" quantum- 
pumping regime. For instance, in mesoscopic structures, 
this regime is attained by modulating electronic poten- 
tials with two low-frequency fields with a relative phasei^ 
More precisely, assuming that the electrons propagate co- 
herently through such a device, it can be shown that heat 
is pumped at a rate proportional to fcsTr^o, while the 
rate at which energy is dissipated scales as TiflQ, where 
T is the temperature of the cold reservoir and flo is the 
frequency of the ac field. As in the cases of charge^ 
and spin^ pumping, coherence is expected to be an im- 
portant ingredient for this type of heat pumping. These 
mechanisms act on the electronic degrees of freedom, and 
therefore they are ineffective in cooling insulators. Cool- 
ing procedures based on phonon manipulation, however, 
apply equally well to metals and insulators. Phononic 
refrigeration is the topic of this paper. 

The minimal microscopic model for acoustic phonons 
consists of masses coupled by springs. Heat transport 
and thermalization in chains of coupled oscillators have 



a rather long history in statistical physics since the pi- 
oneer work by Fermi, Pasta, and Ulam^- Several stud- 
ies of different classical and quantum models have an- 
alyzed the local temperature profile of chains of cou- 
pled oscillators in contact with two reservoirs at differ- 
ent temperatureSfiS, The striking feature found in these 
studies was the violation of the Fourier law, according 
to which the local temperature is expected to drop lin- 
early along the chain. This behavior is the phononic 
counterpart of the abrupt drop of the local voltage 
at the connections between a finite size electronic sys- 
tem to two reservoirs at different chemical potentials, 
which has been characterized by the concept of "con- 
tact resistance" »ii More recently, heat pumping has been 
analyzed through a two-level system driven by an ac 
field coupled to phononic baths at different tempera- 
tures,— between semi-infinite harmonic chains with a 
time-dependent coupling, ^'^ and in anharmonic molecular 
junctions in contact to phononic baths. ^"^ 

In a recent work^i^ we have shown that when a barrier 
(or perturbation) moves along a cavity (which we refer as 
pipe) connecting two gases of acoustic phonons, a sizable 
amount of heat can be transferred from the colder to 
the hotter gas. Assuming that the barrier is perfectly 
rcfiective and its motion is slow enough to allow for a 
fast thermalization of the phonon gas, we can imagine a 
cycle where: (i) a barrier is inserted at the end of the 
pipe close to the cold reservoir; (ii) the barrier is then 
adiabatically translated through the pipe until it reaches 
the end in contact to the hot reservoir; (iii) a second 
barrier is then inserted again near the cold reservoir; and 
(iv) the first barrier is finally removed. The process can 
then be repeated as a cycle. The result is a refrigerator 
which docs not rely on quantum coherence. 

The phonon pump proposed in Ref. |l5| used a ther- 
modynamical formulation for the phonon gas. Here we 
present instead a microscopic model of acoustic phonons 
that implements this cooling scheme. The lattice model 
we consider consists of a series of atoms or molecules with 
identical masses coupled by springs. The moving barrier 
in this microscopic model is realized by modulating in 
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FIG. 1; The one-dimensional microscopic model for acoustic 
phonons consists of a chain of particles (atoms of molecules) 
with equal mass mc coupled by springs with elastic constants 
ko. The semi-infinite left portion of the chain is kept at a 
temperature Tl which is lower than the temperature Tr of 
the semi-infinite right portion. A traveling harmonic potential 
is applied locally to the central portion of the chain. This 
potential restricts the longitudinal motion of the particles - 
see last term on the right-hand side of Eq. ((2|. To perform 
a cooling cycle, the local potential is turned on and moved 
from near the cold to the hot reservoir at a constant speed v. 
When it reaches the hot reservoir, the potential is turned off 
and moved back to near the cold reservoir. 



term provides the contact between the central chain c 
and the L, R reservoirs. 

The masses in the central chain are coupled to their 
neighbors via time-independent springs. An external po- 
tential is applied to the central part of the chain in or- 
der to restrict (i.e., pin) the longitudinal motion of the 
masses. We use a local harmonic approximation for this 
pinning potential. (We defer to Sec. |V] a discussion 
about the possible experimental realizations.) The en- 
suing Hamiltonian reads 
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To illustrate the heat pumping mechanism, we break 
down the spring constant of the local pinning potential 
into a train of pulses with a Gaussian shape. 
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time and space a local pinning potential, which restricts 
the longitudinal motion of the particles. The setup is 
shown schematically in Fig. [1] 

The paper is organized as follows. In Sec. HI] we present 
the microscopic model. In Sec. IIII A[ we show the result 
of numerical computations of the heat currents and use it 
to evaluate the lowest temperature the phonon pump can 
cool down. This result is then compared to the prediction 
made in Ref. [T^ based on a continuum model. In Sec. ITVl 
we provide a detailed description of the non-equilibrium 
Green function methodology used to carry out the com- 
putations. Our conclusions are presented in Sec.|Vl Some 
technical details are presented in the Appendix. 



II. MICROSCOPIC MODEL 

We consider the one-dimensional (Id) system sketched 
in Fig. [1] with a central finite chain of N atoms or 
molecules (which we will refer as "masses" hereafter) 
with identical masses TOc coupled by springs of constant 
kc- This system is connected at its left (right) end to a 
semi-infinite chain L {R) of masses niL ("ifl) coupled by 
springs of constant k^ (kn). These left and right chains 
play the role of reservoirs, which we assume being kept at 
temperatures Tl and Tr, respectively. We only consider 
longitudinal vibrational modes. The ensuing Hamilto- 
nian reads 



H = HL + HR + H,{t)+H,. 



(1) 



The first two terms Hlh represent the reservoirs; the 
third Hc{t) describes the central chain, which serves as 
the pump and has an explicit time dependence; the last 



with fc;*"V) = fci e^["'^''(*""^)l'/'"', TIT <t < {n + 1)t, 
where a is the lattice constant. The barrier speed is v, its 
width is a, and the time it takes to traverse the central 
chain is r = aN/v. It is useful to define two characteristic 
frequencies in the problem: the frequency ujc = y^kQ/rric, 
and the pumping frequency flo — 2t:/t. (We shall work 
on units such that h — ks — ^■) 

The contact between the central system and the reser- 
voirs is described by the Hamiltonian 
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Notice that the central region couples to the left reservoir 
at the site / = 1 and to the right one ai I = N (for later 
use, it is convenient to define l^ = I, In = N). 
The reservoir Hamiltonians are given by 
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with a — L, R, in the limit of Na — ?> cxd (which we shall 
take in due time) . The lattice sites are labeled from right 
to left for the L reservoir and from left to right for the 
R reservoir. It is convenient to express the degrees of 
freedom of the L and R reservoirs in terms of normal 
modes. For open boundary conditions, this corresponds 
to performing the following transformation. 
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where 



n = 0, 
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The corresponding Hamiltonians transform into 
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A. Energy balance 

Following a procedure similar to that used in Ref. Is, 
we consider the variation in time of the energy stored in 
one of the reservoirs. We find 



^-^-J2{t)^patl 



(10) 



where the first term is the heat current flowing from the 
connecting site la of the central chain into the a = L, R 
reservoir, which reads 
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with kn a ^ ko 



j^-^^sm{q2), while 
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is the power due to the time-dependent forces acting at 
the connecting site /„. 

Conservation of energy implies that the total power in- 
vested by the external fields is dissipated into the reser- 
voirs at a rate 
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is the power developed by the external fields at the site 
I averaged over one period, while 
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is the dc component of the heat current at the reservoir 
a. 

Notice that a net amount of work has to be invested 
in order to pump heat from one reservoir to the other. 
Carrying out an analysis similar to that developed in Ref. 
y leads to the conclusion that the rate at which the total 
work done by the external fields is dissipated as heat 
flowing into the reservoirs is proportional to 57q. 



In the case of heat transport induced by a temperature 
gradient between the reservoirs, and in the absence of 
time-dependent perturbations, the heat current ([T5|) can 
be in general expressed as 
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where T{^) is the thermal transmission function of the 
structure, while na{uj) is the Bose-Einstein distribution 
function, which depends of the temperature of the reser- 
voir a. 

We evaluate these physical quantities by employing the 
non-equilibrium Green function formalism. Before de- 
scribing the technical details of the calculations, we first 
present the results that demonstrate the refrigeration ca- 
pabilities of the phonon pump. 



III. PHONONIC REFRIGERATION 
A. Mechanism and minimum cooling temperature 

We now show results for the heat currents correspond- 
ing to a perturbation of the form of a train of traveling 
Gaussian pulses described by the modulation shown in 
Eq. (|31). The main goal here is to compare the results 
obtained for our microscopic phononic system to the pre- 
dictions put forward in Ref. llSl , which were derived for a 
continuum model system with a strong, delta-like trav- 
eling perturbation. As mentioned in Sec. Ill Al and in 
Ref. [3, the rate at which heat is produced is propor- 
tional to r^Q. Therefore, we consider a low driving fre- 
quency Qq — 27r/r in order to ensure a low rate of heat 
dissipation into the reservoirs and to achieve a regime 
where heat can be effectively transported from the cold 
to the hot reservoir. The speed of the barrier is tuned to 
V — aN /t, such that a new pulse appears on the left end 
of the central region every time a pulse moves through 
the right end. 

Before showing the results for the dc heat current, it 
is interesting to analyze the effect of a stationary [v = 0) 
pinning potential on the behavior of the heat transmis- 
sion function T(a;) [see Eq. (H^ j. This is illustrated 
in Fig. [2J where 7~(w) is shown for both the unper- 
turbed chain (fci — 0) and the case with a narrow, time- 
independent Gaussian spatial modulation near the center 
of the chain. It is clear from Fig. [5] that the effect of the 
pinning potential is to decrease the transmission proba- 
bility in the low-energy sector of the phonon spectrum. 
Such an effect is an indication that the pinning poten- 
tial, when stationary at a certain site, acts like a barrier, 
blocking heat transport from the hot to the cold reser- 
voir. 

In Ref. UM we have shown that, besides blocking the 
usual heat flowing in the direction of the heat gradient, a 
moving barrier in a phonon gas also enables refrigeration 
under certain conditions. Figure [3] demonstrates that 
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FIG. 2: The transmission function T{iu) defined in Eq. (I42|) 
for a central system of A'^ = 40 atoms without any bar- 
rier (dashed lines) and with a Gaussian barrier with height 
fci = l.Sfco and a = a (solid lines) centered at I = 20. The 
frequency is shown in units of the chain's characteristic fre- 
quency LJc ~ \/ko/mc. 
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FIG. 3: (Color online) The average energy flowing through 
the contacts in a cycle, J^t, as functions of the temperature 
Tl of the left (cold) reservoir (given in units of ljc)- The 
temperature of the right reservoir is kept fixed Tr = 0.04a;c 
and ilo = 1 X 10~^iJc. The central system contains A'^ = 80 
atoms. The barrier has a height fci — 1.5fco, a — a, and moves 
from left to right with a speed v = aNno/{2n). The arrows 
on the horizontal axis indicate two temperature scales, T* 
and Tmin, defined in the main text. 



this is also the case of the microscopic model considered 
in the present work. We consider a finite-size central 
chain with N = 80 atoms and a low pumping frequency 
r^o = 1 ^ 10~^ ujc- We fix the temperature of the hot 
reservoir at Tn = 0.04 Wc and analyze the behavior of the 
heat current at the contacts between the central system 
and each of the reservoirs. In Fig. |3l we can see that, 
for large temperature differences AT — T^ — Tl , the dc 
heat current at the cold reservoir is positive, indicating 
that the flow enters this reservoir, while it is negative at 
the hot one (the flow exits that reservoir). In this situ- 
ation, the direction of the heat current is ruled by heat 



FIG. 4: (Color online) The relative temperature interval {Tr — 
Tinin)/TR as a function of the right reservoir temperatureTij, 
where Tmin is the minimum temperature of the left reservoir 
for which cooling is possible. Systems of sizes A'^ — 40, 80, 
120, and 160 corresponding to black circles, red squares, green 
down triangles, and blue up triangles, respectively, are con- 
sidered (the solid lines are guides to the eye). For each system 
size, a fitting to Eq. p7[) is performed (dashed lines). Other 
parameters are the same as in Fig. [S] 



gradient (hot to cold). As AT decreases, there is a range 
of temperatures of the cold reservoir, T* < Tl < Tmin, 
where the two currents flowing into the reservoirs be- 
come positive. When the temperature of the cold reser- 
voir overcomes Tmin, J l <^ ^^'^ J r > 0- I^i this regime, 
there is a net heat flow from the cold to the hot reservoir, 
which means that the system operates like a refrigerator 
(i.e., heat flows against the temperature gradient). It 
is also interesting to note that when cooling takes place 

(T > Tmin), the absolute value of the currents Jl,r is 
approximately the same, indicating that the dissipation 
of energy due to the action of the external forces is very 
smaU [see Ecl p^ ]. 

In Ref. [T^, we have shown that an estimate for the 
temperature Tmin is given by 
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where Qb — Xyhv/2TT^^ with being c the sound veloc- 
ity of the system and A being the amplitude of a delta- 
function- like moving barrier [see Eq. (4) in Ref. 1S|. 
That calculation was done in a continuous model assum- 
ing that the moving barrier produced negligible heat. It is 
interesting to check if this prediction accurately describes 
cooling in a microscopic model and still holds when there 
is a nonzero heat generation production by the external 
flelds. In order to build a connection between the param- 
eters used in these two distinct formulations, we notice 
that the sound velocity in the microscopic model is given 
by phonon velocity c — ayfco/m^, while the delta- funtion 
amplitude can be related to the substrate spring constant 
through A = {ki/ko) Cujc, where C is a geometrical factor 
related to the pulse shape. 



Figure |4] shows the dependence of the lowest coohng 
temperature Tmin on the hot reservoir temperature T^ 
for systems of several sizes and a common pumping fre- 
quency f^o = 1 X 1Q~^ iOc- Notice that fixing J7o requires 
varying the barrier speed when considering central chains 
of different sizes. In order to have a meaningful compar- 
ison for the different chain sizes, we present the relative 
difference {TB, — T^i^)/Tji. The dashed lines are obtained 
from Eq. p7)) . The geometrical factor C is used as a 
fitting parameter to adjust Eq. (jl7|) to the numerical 
data. The fitted values for C are size dependent. We 
find C - 3.6, 4.1, 4.3, 4.5, for N = 40, 80, 120, 160, respec- 
tively. In all the plots, Eq. (fTT]) provides a description 
of the numerical data for high enough temperatures Tr. 
As the size of the central chain increases, the range of 
temperatures where Eq. p7p fits well the numerical data 
systematically extends to lower values of Tr. The dis- 
crepancy between Eq. (J17l) and the numerical data for 
low values Tr can be understood by noticing that for 
low temperatures the effect of heating due to the exter- 
nal fields becomes relatively large. Indeed, the effect of 
heating is expected to be stronger for the smaller systems, 
which is indeed where the discrepancy seems larger. In 
summary, as the size of the system increases, the behavior 
of the microscopic model approaches that of an infinite 
continuous system and the effect of heating by the exter- 
nal field becomes relatively insignificant. In the large-size 
limit, the results for the microscopic model tend to coin- 
cide with the analytical prediction of Eq. (IT71) and this 
is precisely the observed behavior in Fig. ^. 



B. Cooling efficiency 



The results of Sec. IIII Al indicate that the model de- 
fined in Sec. |TT]can be indeed regarded as a microscopic 
realization of the phonon pump proposed in Ref. Il5l . 
Thus, let us consider the efficiency of the phonon refriger- 
ator in the continuous, phonon gas model, and comparing 
it with that of a Carnot cycle. 

In Ref. 1 151 we showed that the heat extracted from the 
(cold) L reservoir at the end of one cycle and the heat 
injected on the (hot) R reservoir are given by 
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where Vpipc is the volume swept by the moving barrier 
and e^fi and Pl,b. are, respectively, the energy densities 
and the pressures in both sides of the pump. The work 
done in the cycle is given by 
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The pressure and energy of a phonon gas follow the rela- 
tion p = — {dF/dV)rp = 76, where 7 is the Griineisen 



parameter of the lattice and d denotes the spatial 
dimension~2£. 

The cooling coefficient of performance of the refriger- 
ator is given by 
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We can relate the pressures to the temperatures in the 
phonon gas in d dimensions and one finds (see Refj^) 
that pl/pr = {Tl/Tr)'^+\ Wc therefore arrive at the 
simple form 
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This cooling coefficient of performance is to be compared 
to that of a Carnot cycle, COPcarnot = i-it It ) ■ ^^ ^^ 
useful to express the ratio of these two coefficients of 
performance as 
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COPcarnot 1 + {Tl/Tr) + • • • + [Tl/Tr) 

< 1, for < Tl/Tr < 1 . 

It is thus clear that the refrigerator is generically less 
efficient than the Carnot one. When Tl is close to Tr, 
the efficiency of the phonon refrigerator approaches the 
Carnot limit. But when Tl ^ Tr, and the efficiency 
is rather low compared to that of a Carnot cycle. When 
Tl/Tr = 1/2, COP/COPcarnot=2/3 for one-dimensional 
pipes, and COP/COPcarnot=4/15 for three-dimensional 
systems. 



IV. METHODOLOGY: GREEN FUNCTIONS 
AND DYSON EQUATIONS 

In order to develop an analytical framework where the 
nonequilibrium heat currents in the microscopy phononic 
model can be evaluated, we follow a procedure analogous 
to that introduced in Ref. |l6| for fermionic systems. We 
begin by defining the following lesser and greater Green 
functions, 
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and 
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as well as the retarded one, 

D^AtX) = -tQ{t-t'){[xi{t),xAt')]) 

= e{t~t')[D<,it,t')-D>,it,t% (24) 

For coordinates along the central chain (not including the 
contacts), the Dyson equations for the retarded functions 
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By writing the Dyson equation for D^ along the contacts, 
it is possible to integrate out the degrees of freedom of 
the reservoirs. The result is 
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Here, ry > is an infinitesimal and £"«,« = 
^fcQ(l -cosg^)/ma. 

On the other hand, the Dyson equation for the lesser 
Green functions read 
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where Eq. (|3T|) corresponds to coordinates along the cen- 
tral chain while Eq. p2p corresponds to one coordinate 
along the chain and the other on the reservoir. Here we 
have introduced 
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with na{uj) = l/(e'^/-^° — 1) being the Bose-Einstein dis- 
tribution, which depends on the temperature T^ of the 
reservoir a. 

The strategy now is to solve Eq. (pS)) to evaluate the 
retarded function and then use the result in Eqs. ((3T|) 
and (l32]l to calculate the lesser functions. 



A. Exact solution of the Dyson equation 

The exact functions D^ are obtained by solving the 
linear set of equations represented by Eq. (pS)) . It is also 
convenient to use the representation 
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b"{t,uj) = 7)(")(w)+ ^ e~'^^'''b^\t,uj + k^o) 

X M^')^(o)(tj), (37) 

where we are considering the Fourier representation of 
the time-dependent part of the Hamiltonian, M'^^^(t) = 
^^fcM^^^e-**^""*, while 

£>(•') (w) = [cj2i-Af(")-S^(cj)]-i (38) 

is the retarded Green function corresponding to the cen- 
tral chain of harmonic oscillators connected to the reser- 
voirs but free from the time-dependent modulation of 
the spring constants. In our case, the frequency of the 
ac modulation is Hq = 2tt/t = 2TTNa/v. The exact so- 
lution for this set of coupled linear equations can be ob- 
tained by following the procedure introduced in Ref. Il6| 
for fermionic systems driven by ac potentials. An alter- 
native approach, appropriate for the case of a low fre- 
quency D,Q considered in the present work, is discussed in 
Sec. HYCl 



B. dc heat currents 



1st order 



We now turn to the evaluation of the dc component of 
the heat current flowing from a given reservoir. To this 
end, we first express the heat current in terms of Green 
functions. From Eq. ([TT|) . the dc component of the heat 
flow from the reservoir a can be expressed as foUows, 
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Using Eq. ([32|) and the representation shown in Eq. (p6| , 
we get 
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A solution exact up to 0{VIq) can be obtained by ex- 
panding Eq. (P7)) as follows: 
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,a„i«,„)^i>'»'(-). 

We define the frozen Green function 

Df{t,uj)= [£)(o)(tj)-i-M(i)(i)l~\ 



(43) 



(44) 



in terms of which the exact solution of the Dyson equa- 
tion at 0{V,o) reads 

D(i) {t, Lj) ^ Df {t, Lu) + td^Df {t, uj)^^^^^D<^"^ (uj) . 

(45) 



2nd order 



The first term can be recast using Eq. (IA3|) . leading to 
the following equation for the dc heat current, 






E E 

l3=R.,Lk=-oo 



+ 00 



duj 

2^ 



{uj + kno)[ni3{uj) 



na{uj + kno)]Ta{uj + kno)T0{uj)\Vi^j {k,uj)\ 



(41) 



This equation indicates that a net heat current may exist 
even in the absence of a temperature difference between 
the two reservoirs. Such a fiow contains in general an 
incoming component which accounts for the work done 
by the ac fields and which is dissipated into the reser- 
voirs. As we discuss earlier, under certain conditions, 
a net current between the reservoirs can be established; 
this current can go against a temperature gradient, there- 
fore allowing for cooling. 

In the case of a temperature difference AT between 
the reservoirs, Tr — Tj^ + AT, and in stationary condi- 
tions, i.e. in the absence of ac fields, Eq. (PT|) reduces to 
Eq. (mi) where /°'> = -jjj°' = j(o), being the thermal 
transmission function 



T{u>) = Ta{uj)Tj3{uj) 



D 



(0) 



(c) 



(42) 



C. Low driving frequency solution of Dyson 
equation 

Let us consider the Dyson equation in Eq. (|37l) in the 
limit of low driving frequency flo. 



To obtain the solution exact up to ©(J^q), we consider 
the following expansion of equation Eq. ([37]): 

D{t,uj) - L)(°)(a;)+L)(t,cj)M(i)(i)L)(°)(a;) 

at 
- -dlD{t,uj) ^^/^ J^W(c.). (46) 

The solution exact up to 0{Vt^) is 

m{tM = Df{uJ) + ^^^m\t,u)'^^^DJ(u:) 

^-dlDf{t,oj)^^^Df{oj). (47) 

One can obtain this Green function numerically by first 
discretizing the time in the interval < t < r, then solv- 
ing Eq. (|T7|) for each one time, and finally evaluating the 
Fourier transform to obtain the Floquet representation 
with which the dc component of the heat current can be 
calculated from Eq. (|iT|) . 



D. Low driving frequency expansion of the dc heat 
current 



Let us now use the low-frequency expansion in the ex- 
pression of dc heat current flowing into the reservoir a 
given by Eq. (|4T|) . We first expand na{u} + kflo) and 
Taiijj) in powers of flo- By keeping terms up to 0{D.'^), 
we can identify the following decomposition in the net 
current flowing into the reservoir 






■th 



X+iT'- 



J 



(48) 
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where "thermal" component J^ is due to the gradi- 
ent of temperature and flows from the hot to the cold 

— p 
reservoir, the "pumped" component J^ is purely induced 

by the driving and accounts for the cooling mechanism, 
the "mixed" component J^ results from an interfer- 
ence process of the temperature gradient and the cooling 

mechanism, while the "dissipative" term J^ accounts for 
the dissipation of energy into the reservoir due to the 
action of the ac fields. 

The explicit expression for the thermal component is 



-r=dis 



—th 



— [n^(w) - n„(w)]wr/(w), 

ZTT 



Tfiu^) = r4a.)|P/,,„,,,(o,a;)|2rM^), 



(49) 



with P ^ a, while 2?j^;^.;^(0,aj) is the dc component 
matrix element of the frozen Green function (|44p . Notice 

that J j^ = ~Jji ■ Also, recalling that the left reservoir 

is the cold one, J^ > 0, indicating that this component 
flows from right to left. It is also interesting to note that 
the frozen transmission function corresponds to the chain 
with the barrier. Thus, is practically vanishing for small 
|ci;|, as shown in Fig. [2l 

In order to write an expression for the pumped com- 
ponent, we take as a reference temperature the one cor- 
responding to the cold reservoir, T^, resulting 



j: = 



E 

3=L.R k 



duj 

L 

2-K 



E / ^^^^0 ( - 



did 



Xra,{Lu)\Vf^l^j^{k,i0)\^Tp{uj). 



(50) 



In Appendix [B] we show that J^ = ~Jb.- ^^ '^^'" case, the 
driving mechanism is designed to cool the left reservoir, 
thus, J^ < and this current flows against the temper- 
ature gradient. This component is oc ilo and we neglect 
eventual contributions oc J7q or terms of higher order in 

rio. 

The mixed components read 



-jth~p 
•J L 



7, ^ 






uj) — nL{uj)]ujkno X 



X'/,i,Ar(fc,iu)prfl(tj), 



-T-th—p —th 



-jth-p 
•^ R,0 



= J 



R,0 



-jth-p 
•J J 1 



T D 1.-. ^ 



=L.R k 



dujdT 



AT X 



TLioj)\Vf,,,^{k,w)\'T^{uj) 



(51) 



which are oc JIqAT, being AT = Tj^ — T^. Approxi- 
mating the integrand of the last term, d^nL/{dojdTL) ~ 
dnL/{TLduj), which is the leading contribution in oj/T^, 

we get 7*^^7 ^ AtTj^/Tl - -ATJI/Tl. 

Finally, in J^ , for which we do not provide the explicit 
expression for sake of simplicity, we collect all the terms 



that are oc fl^ and satisfy J'^" > 0, a — L,R. Thus, 
this contribution describe pure dissipation of heat gener- 
ated by the ac fields that flows into the reservoirs. This 
contribution is vanishing small for small Hq. 



SUMMARY AND CONCLUSIONS 



We have considered a microscopic model for the 
phononic refrigeration mechanism proposed in Ref. llSl . 
The model describes a one-dimensional lattice of atoms 
or molecules with identical masses connected by springs 
and in contact with cold and hot reservoirs at its ends. 
Considering only longitudinal vibrational modes, a mov- 
ing "barrier" is realized by applying a local harmonic 
time-dependent pinning potential to the masses. The 
cooling cycle consists of propagating the barrier along the 
chain, going from the cold to the hot reservoir with con- 
stant speed. When it reaches the hot reservoir, the bar- 
rier is returned to near the cold reservoir and the cycle is 
repeated. We have shown that such a perturbation, when 
static, suppresses transmission within the low-energy sec- 
tor of the phononic spectrum. In this sense, our model 
describes a situation similar to that found in electrostric- 
tive polymer systems discussed in Ref. [13, where an elec- 
tric field is used to open a local phonon gap. 

The Hamiltonian described in Eq. ^ could be real- 
ized experimentally in several ways. For instance, in the 
case of a chain consisting of polarizable molecules, the lo- 
cal pinning of intermolecular vibrations could be achieved 
through the application of a inhomogeneous electric field, 
which would be swept along the chain to produce a mov- 
ing phonon barrier. 

We note that the static version of the model described 
in Sec. |ll] is rather standard in the study of heat trans- 
port along harmonic chains with a single mass per unit 
cell, including terms beyond the usual intermass coupling 
but within the harmonic approximation^^ In our case, we 
have considered a time-dependent version of those exten- 
sions in order to model the moving barrier. 

The evaluation of heat currents requires the use of 
a formulation capable of handling both a temperature 
gradient and the time-dependent perturbation simulta- 
neously. To that end, we employed an analytical ap- 
proach based on non-equilibrium Green functions. This 
treatment is similar to the one introduced to study elec- 
tronic systems driven out-of-equilibrium by harmonic 
time-dependent voltages; here, it was suitably adapted 
to phonon propagators j^^ We have also introduced a par- 
ticular strategy that allows for solving Dyson equations 
in the limit of weak driving frequencies and large ampli- 
tudes. The propagators can then be evaluated numeri- 
cally through a time discretization and Fourier transform. 
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Following the same steps of Ref. Isl, we start from 






duji Di^i^{ki,u:i)T a{uJi)D*i, , (fc2,a;i)A<'>(w) 



a; + 177 - (cji + k2^o) 



with A^(a;) — —na{uj) and A^(a;) = — [1 — na{uj)]. Therefore, 

fci a=R,L' 

which implies the following identity: 



A,,(fc,a.)=E E f 



27r 



U! + ij] ~ (wi + fefio) 



fc' a=R,L 



(Al) 



(A2) 



V^'{k,uj)-Vli{-k,u} + kflo)^-iJ2^ E 'Du^ik + k',uj-k'no)Taiuj-k'no)'Dl,iJk',uj-k'nQ). (A3) 



Appendix B: Conservation of the pumped 
component of the net heat current 



and recast Eq. (ISUl) as 



The goal of this appendix is to show that the pumped 
dc heat current ([501 is conserved, in the sense that it 
satisfies 



E^.-o- 



(Bl) 



a=L,R 



T 



T Jo J 2Tr \ duj 



Im 



sit,Lj)dtSHt,Lu) 



(B3) 



To this end, we use the relation between the frozen which has been proved to satisfy Eq. ((BH on the basis of 

the Birman-Krein relation, dln{detS) — Tr(S'(iS'^) (where 
det(X) and Tr(A') are the determinant and the trace of 
a matrix X, respectively), applied to the frozen matrix, 



Green's funcion and the frozen scattering matrix,— 



(B2) which is unitary 
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